Three-body monopole corrections to the realistic interactions 
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It is shown that a very simple three-body monopole term can solve practically all the spectroscopic 
problems — in the p, sd and pf shells — that were hitherto assumed to need drastic revisions of the 
realistic potentials. 

PACS numbers: 21.60.Cs, 21.30.+y, 21.10.-k 
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The first exact Green's function Monte Carlo (GFMC) 
solutions for A > 4 nuclei confirmed that two-body (2b) 
interactions fell short of perfectly reproducing experi- 
mental data [Q. The inclusion of a three-body (3b) 
force lead to excellent spectroscopy, but some problems 
remained for the binding and symmetry energies and 
spin orbit splittings. No core shell model calculations 
(NCSM) Q have recently developped to the point of 
approximating the exact solutions with sufficient accu- 
racy to provide a very important — though apparently 
negative — result in 10 B Q: While in the lighter systems 
the spectra given by a strict two-body potential are not 
always good — but always acceptable — in 10 B, the spec- 
trum is simply very bad. 

My purpose is to show the striking analogy between 
this situation and what occurs in conventional (Ofku) 
shell model calculations with realistic G-matrices, and 
then explain how a very simple 3b term can solve prac- 
tically all the spectroscopic problems — in the p, sd and 
pf shells — that were hitherto assumed to need drastic 
revisions of the realistic (R) 2b potentials. 

The first realistic matrix elements Q, and the first 
large scale shell model codes || appeared almost simul- 
taneously. Calculations for up to five particles in the sd 
shell gave very satisfactory results, but the spectrum of 
22 Na (i.e., (sd) 6 T = 0) was very bad @. (Note that 
10 B is (pfT = 0). At the time nobody thought of 3b 
forces, and naturally the blame was put on the 2b ma- 
trix elements (V s J t ^ v , stuv are subshells). The proposed 
phenomenological cures amounted to fit them to the ex- 
perimental levels. Two "schools" emerged: One proposed 
to fit them all (63 in the sd shell) , and lead eventually to 
the famous USD interaction || [I. The alternative was 
to fit only the centroids, given in Eqs. (0j2|). 
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They are associated to the 2b quadratics in number (n s ) 
and isospin operators (T s ), Eqs. (@,H), and they define 
the monopole Hamiltonian, Eq. (|5j), in which we have 
added the single particle (lb) term. The idea originated 
in Ref. 10 , where it was found that the Kuo Brown 
(KB) interaction in the pf shell Jll| could yield excellent 
spectroscopy through the modifications (/ = fj/2, r = 
J5/2, P3/2, P1/2), 



V£(KB1) = Vf r (KB) - (-) T 300 keV, 
Vff (KB1) = V$ f (KB) - 350 keV, 
Vff (KB1) = V^(KB) - 110 keV. 
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The validity of this prescription was checked in pertur- 
bative calculations |l2| , and convincingly confirmed for 
A — 47 — 52 once exact diagonalizations became feasi- 
ble H 

In what follows / will stand generically for 
(P3/2, d 5 / 2 , /r/2) in the (p, sd, pf) shells respectively. 
Obviously r = P\i% and r = d z / 2 , s i/2 f° r the p and sd 
shells. 

Nowadays the 2b NN potentials are nearly perfect, 
and the calculations are exact. Therefore, the blame for 
bad spectroscopy must be put on the absence of 3b terms. 
Which means that the monopole corrections must be 3b 
and Eq. <M) must be supplemented by 
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where n s t u = n rs n t , or n r (n r — l)(n r — 2)/6 and similar 
forms for T stu . To simplify matters we — tentatively — 
allow only contributions of the type n s t(n — 2) and 
T s t(n — 2), i.e., 2b terms modulated by the total num- 
ber of particles n. It should be borne in mind that a 3b 
interaction also produces 2b pieces in the model space, 
exactly in the same way that the 2b interaction produces 
the single particle splittings by summing over the core 
orbits c of degeneracy D c , 
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Note that a 3b potential will produce both lb and 2b 
terms. We need not worry about the former because 
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FIG. 1: Excitation energies for B referred to the J 
lowest state. See text 
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FIG. 2: Excitation energies for 
lowest state. See text. 
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they correct e s which will be taken from experiment as 
traditionally done. The latter, together with the 3b part 
will transform the realistic (R) 2b centroid V^(R) into 
KI(R1) = V£(R) + (al + Pin) = V s T t (R) + X 'J- 

H m can be characterized by demanding correct sin- 
gle particle and single hole spectra around closed shell 
nuclei 18 . This set (cs ± 1) is taken to include the dif- 
ferences in binding energies (gaps) 2BE(cs) — BE(cs + 
1) — BE(cs — 1). The major monopole correction involves 
the gaps around 12 C, 28 Si, 48 Ca and 56 Ni which are too 
small to produce the observed double magicity pl|. It 
will be taken care by a single linear form n = n(n). The 
generalization of Eq. (0) is then 
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The single particle splittings above the / closures are 
quite well given by some R interactions. Hence the 
corrective term Xf r — which will prove useful in the sd 
shell — is most likely to have a 2b origin, %// ^ s intro- 
duced only for completeness and will be altogether dis- 
regarded. Xrr' mus t play an important role because the 
single hole states (at A — 15, 39 and 79) p3] are severely 
missed. However, they have little influence on the nuclei 
we shall study (at the beginning of the shells) . 

For the 10 B spectrum in Fig. |l| the black squares show 
the results of Navratil and Ormand (NO) j| Fig 4, 67jO] 
for the low lying T = states in 10 B. The black circles 
correspond to the bare KLS G-matrices plj |2^ used in 
p3[ . The agreement with experiment (lines) ]24fl is poor 
but the agreement between the calculations is good. This 
is not a joke, but an important remark: NO provides 
the foundation for a conventional G-matrix study. As 
emphasized over the years [|l0[ ^3[ ^5|, the realistic G- 
matrices are very close to one another and will provide 
good spectra once monopole corrected. Absolute energies 



and strength functions are another matter, and much re- 
mains to be learnt from exact and no-core results. 

The open pentagons in Fig. [l] correspond to the classic 
Cohen Kurath fit Q (CK). The open squares and circles 
refer to the KLS interaction with a k — 1.1 correction in 
Eq. (|sj) . The open squares test the influence of the x^ r 
term through a uniform attraction of 1.5 MeV (in CK 
it is about 3 times as large). Conclusion: there is not 
much to choose between the two LKS corrected cases. 
Moreover, they are practically as good as CK except for 
the second J — 3 level. 

There are two reasons not to dwell any longer in the 
p shell. The first is that the aim of this letter is to show 
that the monopole corrections must be 3b, i.e., k must 
be linear in n, which demands examining cases of suffi- 
ciently different n. Unfortunately this is impossible with- 
out bringing in the other possible contributions: For ex- 
ample, Xrr i s no t ver y significant in 10 B (n = 6), but it 
is important in 12 C (n — 8) and crucial in 14 N (n = 10). 
Therefore, there is no way of exploring what a single term 
in Eq. (||) does: all must contribute. As it happens — and 
this is the second reason — the full exploration has been 
done (2j|, and the results were excellent. At the time, 
the problem was that the 3b contributions turned out to 
be large and important, and the authors did not know 
what to do with them. 

For the 22 Na spectrum in Fig. || the black squares show 
the results for the venerable KB ||. The black circles 
correspond to the BonnC (BC) G-matrices |27], . The 
agreement with experiment (lines) [ p9| is poor but the 
agreement between the calculations is good. Again, this 
is not a joke, but an important remark: as mentioned, 
there are very little differences between the realistic G- 
matrices. The open pentagons correspond to Wilden- 
thal's USD || . The open squares and circles refer to the 
KB and BC interactions with k = 0.9 and 0.85 correc- 
tions respectively. We shall come to the triangles soon. 
Though USD is closer to experiment, the corrected R 
interactions do definitely well. 
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FIG. 3: Excitation energies for 23 Na and 24 Mg referred to the 
J = 3/2 and lowest states repectively. See text. 



The story repeats itself for 23 Na and 24 Mg in Fig. [|. 
The notations are the same as in Fig. | The agreement 
with experiment is now truly satisfactory, and the plot- 
ting technique adopted makes the physics quite evident: 
the trouble with a 2b-only description is that the excited 
band K = 1/2 in 23 Na, and the K = 2 (7) band in 24 Mg 
are too low. 

The open triangles in Figs. || and ||show what happens 
with KB when — instead of keeping k fixed — we increase 
it by steps of 0.5 per n. Though there is an improve- 
ment, it is not sufficient to claim the irrefutability of a 3b 
mechanism. The proof comes when we move to Figs. ^: 
In 27 Si, 28 Si, and 29 Si the local value of k (open squares 
and circles) has decreased to 0.60 for KB and to 0.55 
for BC. A constant n is totally ruled out, while he lin- 
ear law (triangles) does quite well. Clearly, the 3b terms 
are indispensable. The superb 2b-only USD fit was ob- 
tained mostly through the massacre of a strong JT = 20 
pairing term that is a constant feature of the R interac- 
tions, which makes USD R- incompatible 25, Section V]. 
This has been known for some time and it is only occa- 
sionally that trouble may arise. The problem has been 
the difficulty, so far, of obtaining an R-compatible fit of 
comparable quality. The mild exception comes from p3] 
where, as in the p shell, the 3b contributions turned out 
to be so large and important, that the authors did not 
know what to do with them. 
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FIG. 4: Excitation energies for 27 Si, 28 Si and, 29 Si referred to 
the J = 5/2, and 1/2 lowest states respectively. See text. 



In the pf shell KB1 (or KB3) is very good for A = 47- 
52 but it produces too large a gap at 56 Ni (7.5 MeV 
against the observed 6.3 MeV). The most serious prob- 
lem comes from the first BE2(2 ==> 0) transition in 58 Ni 
which falls short of the observed value (140 e 2 fm 4 ) by a 
factor f» 0.4. Here it is expedient to replace the constants 
in Eq. (0) by linear terms that have the same value at 
A — 48, and are reduced by a factor w 0.7 at A = 56 
( Eq. (||) works as well). The situation in 56 Ni becomes 
consistent with experiment but in 58 Ni it remains unac- 
ceptable. The problem is solved by BonnC (BC) [ p7[ , 
with the same k reduction of « 0.7 in going from A = 48 
to A = 56. The key difference between KB anb BC is that 
the intensity of the quadrupole force (in MeV, extracted 
as in pi) is 2.7 for KB and 3.2 for BC. This discrep- 
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FIG. 5: Backbending in 48 Cr. See text. 



ancy is somewhat disturbing, but it does alter the basic 
fact that 3b monopole terms are necessary. Fig. || shows 
that for k = 0.43, BC produces a backbening pattern 
in 48 Cr that is practically as good as the KB3 one. At 
K = 0.28 — the correct value around A = 56 — the agree- 
ment with experiment is destroyed. 

There are several other indications that a 3b interac- 
tion is essential. Perhaps the most significant is the fol- 
lowing: The monopole centroids Vf r ^{sd,) mus t be such 
that when / 7 / 2 fills the d, (I = 2) orbits are depressed 
with respect to the s, (/ = 0) one fL8|| . However, it is 
clear from the spectrum and the spectroscopic factors in 
29 Si that the filling of rf 5 / 2 favours the p, (I = 1) orbit(s) 
over the /, (/ = 3) ones [E9|. A 2b-only assumption leads 
to a contradiction: if / 7 / 2 acting on the sd shell favours 
the larger I orbits, 1^5/2 acting on the pf shell must do 
the same. Without unacceptable ad-hoc assumptions, a 
2b mechanism cannot do otherwise but a 3b one can. 

From what we have seen, 3b monopole forces make 
things simpler, and there are good reasons to believe that 
the formidable task of a full treatment — including multi- 
pole terms — need not be inevitable . A recent generation 
of 3b potentials |50| has made it possible for the exact so- 
lutions to eliminate the more offending aspects of the 2b 
10 B spectrum |3l]| . It will be of much interest to check 
whether the underlying mechanism corresponds to the 
one proposed in this letter. At any rate, a full character- 
ization of the 3b potentials is not an easy matter, and it 
could be hoped that information coming from shell model 
studies may prove valuable. Especially at a time when 
GFMC and no-core calculations have rigorously estab- 
lished the basic reliability of such studies. 

Several obsevations of Alfredo Poves and Frederique 
Nowacki have been of great help. 
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